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The  paper  analyzes  the  damping  of  the  Alfven  surface  waves  via  ion  parallel  viscosity  at  a 
single  magnetic  interface.  The  dispersion  relation  is  obtained  and  characteristics  curves  are 
drawn  for  the  real  and  imaginary  phase  speeds.  The  modes  thus  obtained  are  two  damped 
Alfven  modes  propagating  in  a way  such  that  when  one  mode  dies  new  mode  appears  after 
a certain  propagation  gap.  The  results  are  applicable  for  the  situation  in  solar  wind  at  AU 
for  values  obtained  from  the  spacecraft  data. 


1.  Introduction 

Alfven  waves  have  been  widely  discussed  in  space 
plasma  and  their  ubiquitous  presence  in  the  solar 
wind  was  demonstrated  by  [1].  They  lose  energy 
as  they  propagate  outward  from  the  Sun.  If  we 
consider  solar  wind  as  the  part  of  the  corona  we  can 
undoubtedly  say  that  Alfven  waves  lead  to  coronal 
heating.  The  wave  damping  in  solar  wind  is  not  well 
understood.  In  what  follows  we  study  the  viscous 
damping  of  Alfven  surface  waves  propagating  at  the 
magnetic  interface  in  the  solar  wind.  The  Braginskii 
viscosity  ([2],  [3],  [4]  and  [5])  is  a tensor  with 

terms  770,  771,772,  describing  viscous  dissipation  and 
terms  proportional  to  773  and  774  as  nondissipative 
and  describe  wave  dispersion  related  to  the  finite 
ion  gyroradius.  Since  can  typically  be  of  the 
order  of  106  in  the  solar  wind,  where  u)Ci  is  the  ion 
cyclotron  frequency  and  T{  is  the  ion  collision  time, 
the  term  770  * parallel  viscosity'  becomes  far  more 
important  than  all  other  terms. 


Ti  =B°'Vy' 

(4) 

V • v = 0,  V ■ b = 0, 

(5) 

where  B0  = B0 z;  v and  b are  the  perturbed  veloc- 
ity and  magnetic  field,  and  pT  is  the  total  pressure 
(plasma  and  magnetic). 

Equations  (3)- (5)  can  be  Fourier- analysed,  assuming 
all  perturbed  quantities  a g(x)  exp(~iw>t  4-  ikz)  i.e. 
d/dt  = —iuj^d/dx  ^ 0,8/dy  = 0,  and  djdz  = ikz. 
to  obtain  the  second  order  differential  equation  for 
the  x—  component  of  velocity,  vx , 

d2 

j~2Vx  ~ m2v*  = °>  (6) 

where 
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We  consider  plasma  media  occupying  half  spaces  x < 
0 and  x > 0.  The  solution  to  equation  (6)  in  the 
respective  regions  x < 0 and  x > 0 is  given  by 


2.  Basic  equations 


We  consider  an  incompressible  and  viscous  magneto- 
fluid of  uniform  density  under  the  assumption,  when 
u)cr  1.  Since  the  viscous  dissipation  is  dominated 
by  the  term  proportional  to  ion-parallel  viscosity, 
conventionally  denoted  by  77Q,  the  viscous  term  gets 
reduced  to  ( [6]  and  [7]) 

-77*  V • [(/ - 3BB)BB  : Vv],  (1) 


where  B is  a unit  vector  in  the  direction  of  magnetic 
field.  In  terms  of  cartesian  coordinates  equation  (1) 
simplifies  to 
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(2) 


vXo  — A0em°x,  x < 0, 

(7) 

vxe  = Aee~m‘x,  x > 0. 

(8) 

where  A0  and  Ae  are  arbitrary  constants  and 
subscripts  ’o’  and  ’e1  stand  for  the  regions  x < 0 and 
x > 0 respectively.  We  have  imposed  the  conditions 
mQ  > 0 and  m€  > 0,  and  vx  -y  0 as  x — > ±oc  to 
represent  surface  waves.  Across  the  interface  x = 0, 
we  impose  the  boundary  conditions  that  the  normal 
component  of  velocity  vx  and  total  pressure  pT  (gas 
plus  magnetic)  must  be  continuous. 

3-  Dispersion  relation  and  Discussion 

Using  the  boundary  conditions,  we  derive  the  dis- 
persion relation  in  normalized  form  for  the  surface 
waves 


Taking  the  mass  density  to  be  uniform  and  con- 
stant, the  associated  kinematic  viscosity  coefficient 

iS  Vion  = 1 ]l/p. 

The  relevant  linearized  incompressible  MHD  set  of 
equations  when  the  viscous  dissipation  is  dominative 
by  the  term  proportional  to  ion  parallel  viscosity,  are 

^7  = -VpT  + -p— (B0  • V)b  + Vi0nDvz,  (3) 


(x2  — 1 + 2ixV{)y  r2mQ  — rri2(x2  — aAe  + 2 ixVey , 

(9) 

where 
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Here  ^on,o  and  nt0n,e  arc  the  kinematic  viscosity  co- 
efficients. and  vac  and  va0  are  the  Alfven  velocities 
on  either  side  of  the  interface.  In  the  absence  of  ion- 
parallel  viscosity,  the  dispersion  relation  reduces  to 


_ Po^Ao±PcV%c  nnx 

k2  ~ Po  + p<  U j 

which  is  the  well  known  dispersion  relation  for  Alfven 
surface  waves  in  an  incompressible  fluid.  In  order  to 
know  the  nature  of  waves  in  our  case  we  need  to 
study  the  equation  (9).  On  setting  Ve  = aVG„  where 
a = vJon'c , we  solve  numerically  the  dispersion  rela- 
tion (9)  for  phase  speeds  as  a function  of  Vq  in  the 
context  of  a situation  in  solar  wind  at  1 AU.  The  ra- 
tio of  current  plasma  densities  on  either  side  of  the 
interface.  pofpc  is  taken  0.2.  The  values  of  real  and 
imaginary  phase  speeds  are  obtained  from  the  values 
of  x in  units  of  vao-  Figs.  1 and  2 show  the  variations 
of  real  and  imaginary  phase  speeds  with  the  parame- 
ter Vo  for  r = 0.04,  = 0.44  and  a = 0.5.  It  is  ev- 

ident from  the  figures  that  the  Alfven  surface  waves 
propagating  along  the  interface  are  damped  waves.  It 
is  also  found  that  there  are  only  two  damped  modes 
of  Alfv&i  surface  waves  which  do  not  propagate  si- 
multaneously; the  later  mode  propagates  with  slower 
speed  than  previous  one.  It  is  seen  from  figure  1 
that  the  speed  of  Alfvdn  surface  wave  decreases  as 
ion-parallel  viscosity  increases.  This  wave  becomes 
evanescent  at  a critical  value  of  Vo =0.55.  After  the 
disappearence  of  this  mode  a new  second  mode  arises 
from  0.6  whose  phase  speed  decreases  with  the  in- 
crease in  the  value  of  the  parameter  Vo*  It  is  also 
notable  that  there  is  a small  region  from  Vo  = 0.55 
to  0.6  where  there  is  no  propagation  of  either  of  the 
wave.  This  region  is  called  a non  propagation  region. 
Figure  2 depicts  the  damping  rate  of  the  wave  with 
the  parameter  Vo*  Damping  of  the  mode  increases 
as  Vq  increases  but  decreases  after  the  value  of  0.9. 
Thus  the  modes  of  surface  waves  become  damped 
owing  to  ion-parallel  viscosity  in  an  incompressible 
fluid. 
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Figure  1:  Variation  of  real  phase  speeds  t on/k  with 
parameter  Vq.  Dark  line  represents  first  inode  while 
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